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ABSTRACT 

This paper proposes a class of estimators for population correlation coefficient 
when information about the population mean and population variance of one of the 
variables is not avaliable but information about these parameters of another variable 
(auxiliary) is avaliable, in two phase sampling and analyzes its properties. Optimum 
estimator in the class is identified with its variance formula. The estimators of the class 
involve unknown constants whose optimum values depend on unknown population 
parameters.Following Singh (1982) and Srivastava and Jhajj (1983), it has been shown 
that when these population parameters are replaced by their consistent estimates the 
resulting class of estimators has the same asymptotic variance as that of optimum 
estimator. An empirical study is carried out to demonstrate the performance of the 
constructed estimators. 
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1. Introduction 

Consider a finite population U= {1,2,..,1,..N}. Let y and x be the study and auxiliary 
variables taking values y; and x; respectively for the ith unit. The correlation coefficient 
between y and x is defined by 


Px = Syx /(SySx) (1.1) 
where 
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(xi, yi), 1 = 1,2,...,n; the usual estimator of p,, is the corresponding sample correlation 


coefficient : 
T= Syx /(SxSy) (1.2) 
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The problem of estimating p,, has been earlier taken up by various authors including 
Koop (1970), Gupta et. al. (1978, 79), Wakimoto (1971), Gupta and Singh (1989), Rana 
(1989) and Singh et. al. (1996) in different situations. Srivastava and Jhajj (1986) have 
further considered the problem of estimating p,, in the situations where the information 


on auxiliary variable x for all units in the population is available. In such situations, they 
have suggested a class of estimators for p,, which utilizes the known values of the 
population mean X and the population variance S? of the auxiliary variable x. 

In this paper, using two — phase sampling mechanism, a class of estimators for 
p,, in the presence of the available knowledge (Z and S?) on second auxiliary variable z 


is considered, when the population mean X and population variance S? of the main 
auxiliary variable x are not known. 


2. The Suggested Class of Estimators 
In many situations of practical importance, it may happen that no information is 


available on the population mean Y and population variance S?, we seek to estimate the 
population correlation coefficient p,, from a sample ‘s’ obtained through a two-phase 


selection. Allowing simple random sampling without replacement scheme in each phase, 
the two- phase sampling scheme will be as follows: 


(i) The first phase sample s* {s* cU ) of fixed size n,, is drawn to observe only x in 
order to furnish a good estimates of X and S?. 

(ii) Given s*, the second- phase sample s (scs*) of fixed size n is drawn to 
observe y only. 

Let 

s= Sx, = Wn) Dy. =n) Days? =D), 


We write u =x/x* ,v=s/s” . Whatever be the sample chosen let (u,v) assume values in 
a bounded closed convex subset, R, of the two-dimensional real space containing the 
point (1,1). Let A (u, v) be a function of u and v such that 
h(1,1)=1 (2.1) 

and such that it satisfies the following conditions: 
1. The function A (u,v) is continuous and bounded in R. 
2. The first and second partial derivatives of h(u,v) exist and are continuous and 

bounded in R. 


Now one may consider the class of estimators of p,, defined by 
Pra =r h(u,v) (2.2) 
which is double sampling version of the class of estimators 


T, =r f(u",v’) 


Suggested by Srivastava and Jhajj (1986), where u* =*/X, v* =s?/S? and (x, s2) are 
known. 
Sometimes even if the population mean X and population variance S? of x are 


not known, information on a cheaply ascertainable variable z, closely related to x but 
compared to x remotely related to y, is available on all units of the population. This type 
of situation has been briefly discussed by, among others, chand (1975), Kiregyera (1980 
584). 

Following Chand (1975) one may define a chain ratio- type estimator for p,, as 
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where the population mean Z and population variance S? of second auxiliary variable z 





are known, and 


Zz =(1/m)} z; , s2 =(n, -0'5 b; zy 
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are the sample mean and sample variance of z based on preliminary large sample s of 
size n; (>n). 


The estimator p,, in (2.3) may be generalized as 


zN% IONUL fg N «2 \%4 
x X Sy Z S, 
ôu =r( Z) (= 5 5) (2.4) 


where a,'s (i=1,2,3,4) are suitably chosen constants. 





Many other generalization of ,, is possible. We have, therefore, considered a 
more general class of p,, from which a number of estimators can be generated. 

The proposed generalized estimators for population correlation coefficient p,, is 
defined by 


Pia =r t(u,v,w,a) (2.5) 
where w=2"/Z ,a=s"/S? and t(u,v,w,a) is a function of (u,v, w,a) such that 

t (1,1,1,1)=1 (2.6) 
Satisfying the following conditions: 
(i) Whatever be the samples (s* and s) chosen, let (u,v,w,a) assume values in a closed 
convex subset S, of the four dimensional real space containing the point P=(1,1,1,1). 
(ii) In S, the function t(u,v,w,a) is continuous and bounded. 


(iii) The first and second order partial derivatives of t(u,v,w, a) exist and are 
continuous and bounded in S 
To find the bias and variance of ,, we write 


sy =S)(l+e,),¥=X(+e,),x° =X(1+e;),s2 =S (+e) 
sy =Si(l+e3),7 =Z(+e;),s, =S7(1+e4),5,, =S,,(+e,) 


such that E(e) =E (e))=E(e2)=E(es)=0 and E(e;)=0 Y i= 1,2,3,4, 
and ignoring the finite population correction terms, we write to the first degree of 
approximation 
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non-negative integers. 
To find the expectation and variance of f, we expand ¢t(u,v,w,a) about the point 


P= (1,1,1,1) in a second- order Taylor’s series, express this value and the value of r in 
terms of e’s . Expanding in powers of e’s and retaining terms up to second power, we 
have 


EC Bu =P ton) (2.7) 


which shows that the bias of 6,, is of the order nand so up to order n” , mean square 
error and the variance of fp, are same. 

Expanding (ba =P ie retaining terms up to second power in e’s, taking 
expectation and using the above expected values, we obtain the variance of 6,, to the 
first degree of approximation, as 


Var(6 a) =Var(r) + (p3, /M[C}t (P) + (8o — Dt; (P) — At, (P) — Bt, (P) + 28 o39C yt (Pt (P)] 
-= (pi, [DICE (P) + 6 oo — tz (P) — C2t3 (P) — Boa — 1)t4 (P) - At, (P) - 
Bt, (P) + Dt, (P) + F t, (P) + 28 o39C yt, (P)ta (P) = 28 003C -t3 (P)ta (P)] 
(2.8) 


where ¢;(P), t2(P), t3(P)and t,(P) respectively denote the first partial derivatives of 
t(u,v,w,a) white respect to u,v,w and a respectively at the point P= (1,1,1,1), 


Var(r)= (pix IMEE / p 2, )+ (1/4) oyo +8 499 + 28 999) — {8130 +8310)/ Py] (2.9) 
A= (8 19 +8930 — 26 120 / Pyx )FCx> B = 6 220 +8 040 — 218430 / Py IS 
D= {ò 991 +8091 — 26 114 / Py ECF = 8 02 +8 022 — 26 112 / Py dF 


Any parametric function ¢(u,v,w,a) satisfying (2.6) and the conditions (1) and (2) can 
generate an estimator of the class(2.5). 


The variance of 6,, at (2.6) is minimized for 
[46 040 1) = Bs a] 
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Thus the resulting (minimum) variance of 6,, is given by 
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It is observed from (2.11) that if optimum values of the parameters given by 
(2.10) are used, the variance of the estimator 6,, is always less than that of r as the last 


two terms on the right hand sides of (2.11) are non-negative. 
Two simple functions ¢(u,v,w,a) satisfying the required conditions are 


t(u,v,w,a)= 1+a,(u-1) +a,(v-]l) +a,(w-)+a,(a-]) 
t(u, Vv, W,a) =u"'v? wa 

and for both these functions tı(P) =a,, t2 (P) =a, , ts (P) =a, and t4 (P) =a,. Thus one 
should use optimum values of a,,a,, aand a,in 6,, to get the minimum variance. It is 


to be noted that the estimated 6,, attained the minimum variance only when the optimum 


values of the constants a, (i=1,2,3,4), which are functions of unknown population 
parameters, are known. To use such estimators in practice, one has to use some guessed 
values of population parameters obtained either through past experience or through a 
pilot sample survey. It may be further noted that even if the values of the constants used 
in the estimator are not exactly equal to their optimum values as given by (2.8) but are 
close enough, the resulting estimator will be better than the conventional estimator, as has 
been illustrated by Das and Tripathi (1978, Sec.3). 

If no information on second auxiliary variable z is used, then the estimator f, 
reduces to 9, defined in (2.2). Taking z = 1 in (2.8), we get the variance of 6,, to the 
first degree of approximation, as 
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Thus the minimum variance of §,,, is given by 
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It follows from (2.11) and (2.14) that 
2 D m2 
min. Var( a )-min.Var( 6 ,,)=(p2,/m) [ 2 + OCR FF (2.15) 
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which is always positive. Thus the proposed estimator Ĥ„ is always better than 6,,. 


3. A Wider Class of Estimators 
In this section we consider a class of estimators of p,, wider than (2.5) given by 


Pea =8(1,UV,W, a) (3.1) 


where g(r,u,v,w,a) is a function of 7,u,v, w,a and such that 
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Proceeding as in section 2, it can easily be shown, to the first order of approximation, that 
the minimum variance of f ga is same as that of Oia given in (2.11). 
It is to be noted that the difference-type estimator 
r= r+ a, (u-l) +a, (v-1) + a, (w-1) + a, (a-1), is a particular case of Poa , but it is 


not the member of P in (2.5). 


4. Optimum Values and Their Estimates 


The optimum values ¢;(P) = a, t(P) = B ,t3(P)=y and ¢t,(P) =6 given at 
(2.10) involves unknown population parameters. When these optimum values are 
substituted in (2.5) , it no longer remains an estimator since it involves unknown 
(a, B Yy, ), which are functions of unknown population parameters, say,,5 „ım (P, q,m= 


0,1,2,3,4), Cv C, and p,, itself. Hence it is advisable to replace them by their consistent 


estimates from sample values. Let (à, p.f „Š ) be consistent estimators of t)(P),t2(P), 
t3(P) and t4(P) respectively, where 
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We then replace (a ,B ,y ,5 ) by (d, 8,7, ) in the optimum f, resulting in the estimator 
p4 say, which is defined by 


A A 


by =rt (u,v, w,a,0, B78), (4.2) 


where the function t*(U), U= (u,v, w,a,0,B,7 Š ) is derived from the the function 
t(u,v,w,a) given at (2.5) by replacing the unknown constants involved in it by the 
consistent estimates of optimum values. The condition (2.6) will then imply that 
t*(P*)=1 (4.3) 
where P*=(1,1,1,1,a,B 7 ,6 ) 
We further assume that 


sen- ay, npn- -5 4a 
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Expanding t*(U) about P*= (1,1,1,1, a ,ßB ,y ,6 ), in Taylor’s series, we have 


Bia =e (P) + (u = Dti (P*) +0- Dt (P) + (w- Dts (P ) + (a -D (P) + — ats (P’) 
+ (É — BJE (P*) + (y —y)e (p*)+ (5 —5)t; (P*) + second order terms] 
(4.5) 
Using (4.4) in (4.5) we have 


pi =r[l+(u—Da +(v-lB +(w—Dy +(a—1 + second order terms] (4.6) 


Expressing (4.6) in term of e’s squaring and retaining terms of e’s up to second degree, 
we have 


a * 1 * * k * 
(Pia - Pix) = P wl Qes —€) —€,) +a(e, —e,)+ Ble, —e,) +7 e; +ò eq) (4.7) 


Taking expectation of both sides in (4.7), we get the variance of 6,, to the first degree of 
approximation, as 
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which is same as (2.11), we thus have established the following result. 


(4.8) 
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Result 4.1: If optimum values of constants in (2.10) are replaced by their consistent 
estimators and conditions (4.3) and (4.4) hold good, the resulting estimator 6,, has the 
same variance to the first degree of approximation, as that of optimum f . 


Remark 4.1: It may be easily examined that some special cases: 
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(iii) f% =r[l+a(u-l) + B(u-1) +f (w-1) +8 (a-D] 





(iv) Pua = -&(w 1D) - BQ -D-¥(w-) -8(a-DT" 
of f% Satisfy the conditions (4.3) and (4.4) and attain the variance (4.8). 


Remark 4.2: The efficiencies of the estimators discussed in this paper can be compared 
for fixed cost, following the procedure given in Sukhatme et. al. (1984) and Gupta et. al. ( 
1992-93). 


5. Empirical Study 
To illustrate the performance of various estimators of population 
correlation coefficient, we consider the data given in Murthy [1967, P.226]. The variates 
are: 
y=output, x=Number of Workers, z =Fixed Capital 
N=80, n=10, n;-25, 
X = 283.875, Y =5182.638, Z =1126, C, = 0.9430, C, = 0.3520, C, = 0.7460, 
003 = 1-030, Sy, = 2.8664, 5,., =1.1859, S =3.1522, o =1.295, Sg =3.65, 
io. = 9.7491, ôv = 0.9145, 6,,, = 0.8234, ô z0 = 2.8525, 
to = 2.5454, 8 549 = 0.5475, 8 a9 = 2.3377, 555, = 0.4546, Ö-o = 2.2208, 555, = 0.1301, 
so = 2-2667, p, = 0.9136, pẹ = 0.9859, p,, = 0.9413. 


The percent relative efficiencies (PREs) of Pia, Pras Oa With respect to conventional 
estimator r have been computed and compiled in Table 5.1. 


ô 
ô 
ô 
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Table 5.1: The PRE’s of different estimators of p „ 





Estimator r Pra Pa (OF Pia) 




















PRE(.,r) 100 129.147 305.441 





Table 5.1 clearly shows that the proposed estimator 6,,(or p% ) is more efficient 
than rand p,,. 
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